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Abstract 

NASA has conducted a series of Computa- 
tional Aeroacoustics (CAA) Workshops on Bench- 
mark Problems to develop a set of realistic CAA 
problems that can be used for code validation. In 
the Third (1999) and Fourth (2003) Workshops, the 
single airfoil gust response problem, with real geom- 
etry effects, was included as one of the benchmark 
problems. Respondents were asked to calculate the 
airfoil RMS pressure and far-held acoustic intensity 
for different airfoil geometries and a wide range of 
gust frequencies. This paper presents the validated 
solutions that have been obtained to the benchmark 
problem, and in addition, compares them with clas- 
sical hat plate results. It is seen that airfoil geometry 
has a strong effect on the airfoil unsteady pressure, 
and a significant effect on the far-held acoustic in- 
tensity. Those parts of the benchmark problem that 
have not yet been adequately solved are identified 
and presented as a challenge to the CAA research 
community. 

Introduction 

The validation of increasingly sophisticated 
CAA codes is one of the major challenges facing 
the CAA research community today. Recognizing 
this challenge, NASA has sponsored a series of Com- 
putational Aeroacoustics Workshops on Benchmark 
Problems [1-4], with the express purpose of devel- 
oping a realistic set of validated CAA solutions that 
can be used for code validation. 

A fundamental problem that has been included 
in all the Workshops to date is the gust response 


problem. Initially, the focus was on flat plates, both 
for single airfoils and cascades. More recently, the 
focus has been on airfoils with real geometries. In- 
cluded in this latter category is a problem that was 
submitted at the Third [3] and Fourth [4] Work- 
shops, in which respondents were asked to calculate 
the aeroacoustic response of a 12% thick Joukowski 
airfoil to a periodic gust (for different airfoil geome- 
tries and gust frequencies). 

The purpose of this paper is to present the val- 
idated solutions that have been obtained, and to 
identify those parts of the benchmark problem that 
have not yet been adequately solved. We also com- 
pare the validated results with those of classical thin 
airfoil theory, and show that airfoil geometry has a 
strong effect on the airfoil unsteady pressure, and a 
significant effect on the far-held acoustic intensity. 


Benchmark Problem 

Consider the airfoil configuration shown in Fig- 
ure 1. The airfoil is two-dimensional with chord 
length c and angle of attack a. The upstream veloc- 
ity is 

U = Uoc'i + a cos[k ■ (x — iUoot)} (1) 

where x = ( Xi,X 2 ) denotes the spatial coordi- 
nates, k = (^ 1 ,^ 2 ) is the wave number vector, 
a = ( 01 , 02 ) is the gust amplitude vector with 
ai = — eUoo k. 2 /\k\, 02 = e f/oo fci/|fc|, and e is 

a small parameter satisfying e <C 1. Note that 
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Figure 1 Airfoil in a two-dimensional gust. 

a and k satisfy a ■ k = 0, and that k specifies the 
direction of gust propagation. 

The governing equations are the 2-D Euler 
equations 

dp d d 

di + T X {PU) + dy {pv) = ° (2) 

§j.(P u ) + ^(P^+P) + = 0 ( 3 ) 

J^(H + ^ X (P UV ) + J^(pu 2 +p) = 0 ( 4 ) 

+ ^[(£t +P)u] + Wy [{E t +p)v] = 0 (5) 

where p, it, v, p and E t denote the fluid density, 
x and y velocity, pressure, and total internal en- 
ergy per unit volume. The system is closed using 
the perfect gas equation of state. Flow variables are 
nondimensionalized as follows: 


Xi,x 2 

by 

c 

2 

U = (u,v) 

by 

u x 

Co (sound speed) 

by 

u x 

P 

by 

Poo 

P 

by 

PooUl 

T (temperature) 

by 

Too 

t (time) 

by 

c 

2U OQ 

lo = kiUoc (gust angular 

by 

2£/oo 

c 

frequency) 



ki,k 2 

by 

2 

c 


The nondimensional parameter = k\ is 

called the reduced frequency. 

For the following two cases, respondents are 
asked to solve the gust response problem for a 
Joukowski airfoil in a two-dimensional gust with 
k 2 = ki, for reduced frequencies fci = 0.1, 1.0, 
2.0, and 3.0. The nondimensional upstream ve- 
locity is U = i + ea cos (k ■ x — kit), where 
a = (ai,a 2 ) = (-yr.yr). Take e = 0.02. 

For Case 1, the airfoil has a 12% thickness ratio, 
free stream Mach number = 0.5, angle of attack 
a = 0°, and a camber ratio of zero. 


For Case 2, change a to 2° and the camber ratio 
to 0.02. 

The airfoil geometries can be generated as fol- 


lows. Set 

+ 

o 

II 

(6) 

where 

Co' = -ei + «e 2 

(7) 


is a complex constant. Letting z = x + iy denote 
the airfoil coordinates in the complex 2 plane, the 
transformation 


z=( Ci + £)e- ia (8) 

transforms the (i circle defined by equation (6) into 
the desired airfoil shape. 



Figure 2.a Symmetric Joukowski airfoil with 12% 
thickness ratio and zero degree angle of attack. (Case 1) 



Figure 2.b Joukowski airfoil with 12% thickness ratio, 
2% camber, and 2 degree angle of attack. (Case 2) 


For Case 1, use ro = 0.54632753, ei = 

0.05062004, e 2 = 0, d 2 = 0.24572591, a = 0. All 
parameters are specified to eight decimals to ensure 
that the airfoil geometry is precisely defined. Dis- 
cretize the Ci circle in 0, starting from 0 and going 
to 2 7 r, and then apply equation (8) to get the airfoil 
coordinates. The values 0 = 0 and 9 = 2ir map into 
the trailing edge point. 

For Case 2, use ro = 0.54676443, ei = 

0.05062004, e 2 = 0.02185310, d 2 = 0.24572591, a = 
0.034906585. Discretize the £i circle in 0, starting 
from 0 = — / 3 and going to 0 = 2 7r — (3, where 
f3 = 0.039978687, and then apply equation (8) to 
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get the airfoil coordinates. The values 9 = — (3 and 
9 = 2 7T — (3, map into the trailing edge point. 

The above procedure will generate a Joukowski 
airfoil of chord length 2, situated very nearly be- 
tween x = — 1 and x = 1. The airfoil geometries are 
shown in Figure 2. 

For both Case 1 and Case 2, respondents are 
asked to march the discrete equations in time until 
the solution becomes periodic. On the airfoil sur- 
face, calculate the RMS pressure \J ( p ') 2 . In the far 

field, calculate the intensity ( p ') 2 on a circle centered 
about the airfoil center for the following three circle 
radii: (i) R = 2 (one chord length); (ii) R = 4 (two 
chord lengths); (iii) R = 8 (four chord lengths). 


Linearized Benchmark Solution 

Since the gust amplitude is small, the solution 
can be obtained from the linearized Euler equations. 
Let the flow held be decomposed according to 


U (x, t) = U 0 {x) + u(x, t) 

(9) 

p(x,t) = p 0 (x) + p'(x,t) 

(10) 

p(x,t) = p 0 {x) + p'(x,t) 

(11) 

s(x, t) = Sp + s'{x, t) 

(12) 


where the entropy s 0 is constant, and u, p ' , p' , and 
s' are the unsteady perturbation velocity, pressure, 
density and entropy, respectively. Mean how quanti- 
ties are denoted by zero subscripts and are assumed 
to be known. 

Substituting (9) - (12) into the nonlinear Euler 
equations and neglecting products of small quanti- 
ties, one obtains the linearized equations 


DqP' 

Dt 


+ //V • U 0 + V • ( p 0 u ) 


= 0 


Po( 


D 0 u 

~Dt 


+ u-VU 0 )+p'U 0 -VU 0 


D 0 s' 

Dt 


= 0, 


-Vp' 


(13) 

(14) 

(15) 


where ^ + U 0 ■ V is the convective derivative 

associated with the mean how. 

In the present problem, the mean velocity U 0 
can be expressed as the gradient of a potential 4> 0 . 
Therefore, equations (13) - (15) can be reduced to 
the single convective wave equation [5,6] 


D 0 1 D 0 (j) 
Dt { cp 2 Dt ’ 


— V • (p o V0) 
Po 


1 

Po 


y-(p 0 u {R) ), 

(16) 


where the unsteady velocity is decomposed into a 
known vortical component u ^ and an unknown po- 
tential component V0, 

u(x, t) = + V(^. (17) 

The unsteady pressure is given by 

p = -p 0 {x)^-. (18) 

An unsteady aerodynamic code, called 
GUST3D [7], has been developed to solve equation 
(16) for hows with periodic vortical disturbances. 
The code uses a frequency domain approach with 
second-order central differences and a pressure ra- 
diation condition in the far held [8]. GUST3D re- 
quires as input certain mean how quantities which 
are calculated separately by a potential how solver. 
This solver calculates the mean how for interior mesh 
points using a Gothert’s Rule approximation [7]. On 
the airfoil surface, where the Gothert’s Rule approx- 
imation is no longer valid, GUST3D uses the poten- 
tial how solution calculated by FL036 [9]. 

Recently, GUST3D was upgraded to include a 
domain decomposition approach [10]. This has led 
to significant improvements in accuracy. Both ver- 
sions of the code were used to obtain the benchmark 
solutions presented at the Fourth Workshop and in 
this paper. 

Finally, we employ a Kirchoff method [11] to 
obtain the intensity results on the circle of radius 
four chord lengths. The Kirchoff surface is the cir- 
cle of radius two chord lengths. This approach has 
proven to increase accuracy in far-held pressure cal- 
culations. 


Benchmark Results 

In Figs. 3 - 19 we present the benchmark results. 
Figures 3-9 present the RMS airfoil pressure, and 
Figs. 10 - 19 present the acoustic intensity on circles 
of radius one and four chord lengths. For brevity, 
we omit the intensity results for the circle of radius 
two chord lengths. These results essentially paral- 
leled the other intensity results. The hat plate RMS 
results shown in Figs. 3.b - 9.b were obtained from 
a Possio solver [12,13], and the hat plate intensity 
results shown in Figs. 10. b - 19. b were obtained by 
the method of Atassi et al [13]. In each figure, the 
hat plate solution is shown versus a representative 
benchmark solution. Only one benchmark solution 
is shown for clarity. 

The benchmark results selected for each figure 
represent the best solutions available from the Third 
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and Fourth Workshops. Along with the GUST3D 
code, the results are from the PFC6 code of Hixon et 
al [14-16], the CE/SE code of Wang et al [17,18,4], 
the DSEM code of Rasetarinera et al [19,20], the 
BASS code [21-23,4], and the STMA code [24,4], 
It is worth noting that the above codes are all 
based on fundamentally different algorithms, except 
for BASS, which is an upgraded version of PFC6. 
All the codes, except GUST3D, are nonlinear Euler 
solvers. The agreement between the codes ranges 
from good to excellent for most of the figures. It 
should be noted that when results from a particular 
code are not shown, it means either no solution was 
available or that the solution did not agree well with 
the other codes. 

From Figs. 3.b - 9.b, it is clear that airfoil geom- 
etry has a very strong effect on the unsteady pres- 
sure. The flat plate solutions differ markedly from 
the Joukowski airfoil solutions in every case. For 
the flat plate, the strong leading edge singularity is 
clearly visible, whereas for the thick airfoil, there 
is only a pressure spike due to the rounded leading 
edge. Note that the flat plate RMS pressure is iden- 
tical on both sides of the airfoil, whereas the RMS 
pressure for the Joukowski airfoil differs from one 
side to the other. This is due to the airfoil thickness 
and the fact that the gust propagates at a 45° angle 
i.e., k 2 = k t . 

From Figs. 10. b - 19. b, one observes substan- 
tial differences between the intensity profiles for the 
Joukowski airfoils versus the flat plate. As ex- 
pected, the differences tend to diminish somewhat 
on the four-chord-length circle, where the mean flow 
is nearly uniform. However, even at this distance, 
there are still substantial differences for some cases. 

The tables below summarize the benchmark 
problem and the validated solutions that have been 


obtained. In spite of the rather large number of 
validated solutions, it is clear from Tables I and II 
that more work remains to be done. The high fre- 
quency cases k\ = 2.0 and 3.0 indeed pose a signif- 
icant challenge, for a number of reasons. First, the 
incoming gust must be propagated accurately all the 
way from the upstream boundary to the airfoil. For 
the high frequency/short wavelength cases, this re- 
quires a fine mesh in the far field, leading to many 
more mesh points and much longer solution times. 
Second, radiation boundary conditions typically lose 
accuracy at the higher reduced frequencies, making 
it more difficult to obtain an accurate solution. Fi- 
nally, for high reduced frequency cases, there is a 
stronger interaction between the mean flow and the 
convected gust, leading to a more challenging phys- 
ical problem. 

Conclusion 

This paper presents a new set of validated CAA 
solutions to the single airfoil gust response problem. 
The results were obtained from NASA’s Third and 
Fourth Computational Aeroacoustics Workshops on 
Benchmark Problems held in 1999 and 2003. All to- 
gether, the solutions come from six different CAA 
codes, five of which are full Euler solvers. Good 
to excellent agreement is demonstrated in the air- 
foil RMS pressures and far-field acoustic intensities 
for different reduced frequencies and airfoil geome- 
tries. Comparison with classical flat plate solutions 
shows that airfoil geometry has a strong effect on the 
RMS airfoil pressure, and a significant effect on the 
far-field acoustic intensity. Further work is needed 
to obtain accurate results for the high reduced fre- 
quency cases h = 2.0 and 3.0. These cases are pre- 
sented to the CAA research community as a chal- 
lenge yet to be completed. 
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AF Pressure 

Intensity, R=lc 

Intensity, R—Ac 

ki 




0.1 

GUST3D.PFC6, 
CE/SE, STMA 

GUST3D.BASS, 
CE/SE, STMA 

GUST3D.BASS, 
CE/SE, STMA 

1.0 

GUST3D.DSEM, 

CE/SE 

GUST3D.BASS, 
CE/SE, STMA 

GUST3D.STMA, 
CE/SE, DSEM 

2.0 

GUST3D.BASS 

xxxxxxxxxx 

XXXXXXXXXX 

3.0 

GUST3D.DSEM 

CE/SE 

CE/SE 

DSEM, CE/SE 

Table I 

Solution matrix for Case 1. 




AF Pressure 

Intensity, R=lc 

Intensity, I?=4c 

ki 




0.1 

GUST3D,PFC6, 
CE/SE, STMA 

GUST3D.BASS, 
CE/SE, STMA 

GUST3D,BASS, 
CE/SE, STMA 

1.0 

DSEM, CE/SE 

BASS, CE/SE, 
STMA 

BASS, CE/SE, 
STMA, DSEM 

2.0 

GUST3D.BASS 

XXXXXXXXXX 

XXXXXXXXXX 

3.0 

XXXXXXXXXX 

xxxxxxxxxx 

xxxxxxxxxx 

Table II 

Solution matrix for Case 2. 
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Figure 3.b Symmetric airfoil solution (GUST3D) 
versus flat plate solution, k,=k 2 =0.1 
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Figure 4. a RMS pressure on airfoil surface, 
symmetric airfoil, k 1= k 2 =1.0 
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Figure 4.b Symmetric airfoil solution (GUST3D) 
versus flat plate solution, k 1= k 2 =1 .0 



Figure 5. a RMS pressure on airfoil surface, Figure 5.b Symmetric airfoil solution (GUST3D) 

symmetric airfoil, k 1= k 2 =2.0 versus flat plate solution, k 1= k 2 =2.0 



Figure 6. a RMS pressure on airfoil surface, 
symmetric airfoil, ^=^=3.0 


Figure 6.b Symmetric airfoil solution (DSEM) 
versus flat plate solution, k 1= k 2 =3.0 
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Figure 7. a RMS pressure on airfoil surface, 
cambered airfoil, k 1= k 2 =0.1 



Figure 8. a RMS pressure on airfoil surface, 
cambered airfoil, k 1= k 2 =1 .0 



Figure 9. a RMS pressure on airfoil surface, 
cambered airfoil, k 1 =k z =2.0 
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Figure 7.b Cambered airfoil solution (GUST3D) 
versus flat plate solution, k 1= k 2 =0.1 



L.E. x T.E. 

Figure 8.b Cambered airfoil solution (DSEM) 
versus flat plate solution, k 1= k 2 =1 .0 
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Figure 9.b Cambered airfoil solution (GUST3D) 
versus flat plate solution, k 1= k 2 =2.0 
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Figure lO.a Acoustic intensity on circle R = lc, 
symmetric airfoil, kj=k 2 =0. 1 



|p| A 2 cos(0) 


Figure lO.b Symmetric airfoil solution (GUST3D) 
versus flat plate solution on circle R=lc, k =k^=0.1 
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Figure 1 1 .a Acoustic intensity on circle R=4c, 
symmetric airfoil, kj=k 2 =0.1 



|p| A 2 cos(0) 


Figure 12.a Acoustic intensity on circle R=lc, 
symmetric airfoil, kj=k 2 =1.0 
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Figure 1 l.b Symmetric airfoil solution (GUST3D) 
versus flat plate solution on circle R=4c, kj=k 2 =0.1 
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Figure 12.b Symmetric airfoil solution (BASS) 
versus flat plate solution on circle R=lc, k 1 =k 2 =1.0 
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Figure 13. a Acoustic intensity on circle R=4c, 
symmetric airfoil, ^=^=1.0 


Figure 13.b Symmetric airfoil solution (GUST3D) 
versus flat plate solution on circle R=4C, k 1 =k 2 =1.0 



Figure 14.a Acoustic intensity on circle 
R=lc, symmetric airfoil, k 1 =k 2 =3.0 



Figure 14.b Symmetric airfoil solution (CE/SE) 
versus flat plate solution on circle R=lc, k 1 =k 2 =3.0 




Figure 15. a Acoustic intensity on circle R=4c, 
symmetric airfoil, k J =k 2 =3.0 


Figure 15.b Symmetric airfoil solution (DSEM) 
versus flat plate solution on circle R=4c, k 1 =k 2 =3.0 
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Figure 16. a Acoustic intensity on circle Figure 16.b Cambered airfoil solution (GUST3D) 

R=lc, cambered airfoil, k 1 =k 2 =0.1 versus flat plate solution on circle R=lc, kj=k^=0. 1 




Figure 17. a Acoustic intensity on circle Figure 17.b Cambered airfoil solution (CE/SE) 

R=4c, cambered airfoil, k | =k 2 =0.l versus flat plate solution on circle R=4c, k^k^O. 1 
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Figure 18. a Acoustic intensity on circle 
R=lc, cambered airfoil, k ] =k 1 =1.0 


Figure 18.b Cambered airfoil solution (BASS) 
versus flat plate solution on circle R=lc, k ] =k 2 =1.0 




Figure 19. a Acoustic intensity on circle Figure 19. b Cambered airfoil solution (STMA) 

R=4c, cambered airfoil, k =k.=1.0 versus flat plate solution on circle R=4c, k ( =k ,= 1 .0 
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